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The proposal of quantum lithography [Boto et al, Phys. Rev. Lett. 85, 2733 (2000)] is studied via a rigorous 
formalism. It is shown that, contrary to Boto et a/.'s heuristic claim, the multiphoton absorption rate of a 
"NOON" quantum state is actually lower than that of a classical state with otherwise identical parameters. The 
proof-of-concept experiment of quantum lithography [D'Angelo et al, Phys. Rev. Lett. 87, 013602 (2001)] is 
also analyzed in terms of the proposed formalism, and the experiment is shown to have a reduced multiphoton 
absorption rate in order to emulate quantum lithography accurately. Finally, quantum lithography by the use of a 
jointly Gaussian quantum state of light is investigated to illustrate the trade-off between resolution enhancement 
and multiphoton absorption rate. 
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I. INTRODUCTION 

Optical lithography, the process in which spatial patterns 
are transferred via optical waves to the surface of a substrate, 
has been hugely successful in the fabrication of micro and 
nanoscale structures, such as semiconductor circuits and mi- 
croelectromechanical systems. Conventional lithography can- 
not produce features much smaller than the optical wave- 
length, due to the well known Rayleigh resolution limit [ 1 ] . 
As a result, beating the resolution limit for lithography has be- 
come an important goal in the field of optics, with far-reaching 
impact on other research areas including semiconductor elec- 
tronics and nanotechnology. 

Among the many candidates proposed to supersede con- 
ventional lithography, the use of extreme ultraviolet light for 
lithography [2] has met numerous technical difficulties such 
as optics imperfections and photoresist limitations. Other pro- 
posals involve multiphoton exposure JH 0, & 01, so that 
one can still use the more robust optics for long-wavelength 
light, while obtaining some of the resolution improvements 
associated with higher harmonics. The feature size reduction, 
however, is not nominal without taking care of the residual 
long-wavelength features in a multiphoton absorption pattern. 
Yablonovitch and Vrijen proposed the use of several frequen- 
cies and narrowband two-photon absorption to suppress such 
long-wavelength features [3], while a much more radical pro- 
posal by Boto et al. suggests the use of -photon quantum 
interference of N spatially entangled photons |4|]. The so- 
called quantum lithography has several appeals, such as ar- 
bitrary quantum interference patterns, generalization to arbi- 
trary number of photons, and the promise of multiphoton ab- 
sorption rate improvement via the use of entangled photons. 
Hence, despite practical issues such as difficulties in generat- 
ing a high dosage of the requisite entangled photons and find- 
ing a suitable multiphoton resist, interest in quantum lithogra- 
phy has been significant [8]. 

Crucial to the future success of quantum lithography is 
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the supposed enhancement in the multiphoton absorption 
rate when entangled photons are used. The absorption rate 
improvement for frequency anticorrelated photons has been 
proved by Javanainen and Gould BSD and Perina et al. iflQIl . 
Boto et al. further claimed that the absorption rate should also 
improve for the spatially entangled photons used in quantum 
lithography. This promise is absolutely vital to the practi- 
cality of quantum lithography, because, as Boto et al. men- 
tioned, classical multiphoton lithography is already infeasi- 
ble for large N, and quantum lithography would have been 
even worse if the absorption rate was not enhanced, due to 
the much less efficient generation of nonclassical light. Boto 
et al. supported their claim by arguing heuristically that the 
entangled photons are constrained to arrive at the same place 
and at the same time. This argument with respect to the spa- 
tial domain has, however, not been substantiated with a more 
rigorous proof, and has been subject to criticism [11]. Unfor- 
tunately, the time domain treatment |9|, [l0|] cannot be directly 
carried over to the spatial domain, because the former assumes 
a nearly resonant multiphoton absorption process and does not 
require any temporal resolution, but for quantum lithography 
the material response needs to be spatially local to produce a 
high spatial resolution. 

On the other hand, Boto et al. 's formalism contains several 
crucial approximations that remain to be justified. First, the 
photons are implicitly assumed to arrive from a monochro- 
matic source with a well defined free-space wavelength X, 
but the usual quantization method of optical fields consid- 
ers photons as quanta of electromagnetic mode excitations, in 
which frequency appears only as a dependent variable of the 
wave vector. It remains a question how monochromatic opti- 
cal fields as a boundary condition should be treated in quan- 
tum optics. Second, they treat the two optical beams with 
opposite transverse momenta as two discrete modes of optical 
fields, but in free space, transverse momentum is a continu- 
ous variable, so the discrete modes are evidently an approx- 
imation. Third, when discussing the multiphoton absorption 
rate, Boto et al. regards photons as objects in space and time 
that probabilistically arrive at the photoresist, although it is 
well known that this interpretation of photons is fundamen- 
tally flawed 13. 
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In this paper, starting from basic principles, I shall first ex- 
plicitly quantize the electromagnetic fields in Sec.[TTl using as- 
sumptions consistent with Boto et al.'s proposal. The formal- 
ism rigorously treats approximately monochromatic optical 
fields as a boundary condition in the continuous Fock space, 
and hence provides a theoretical underpinning to the proposal 
of quantum lithography. The formalism also shows that, re- 
gardless of the nonclassical spatial properties of the photons, 
there exists an upper bound of multiphoton absorption rate, 
which rules out any significant enhancement of multiphoton 
absorption rate due to spatial effects only. Next, using the de- 
veloped formalism, I shall analyze in Sec.|III]the multiphoton 
absorption rate of the so-called NOON state, and compare it 
with that of a classical state with otherwise identical param- 
eters. The analysis shows that, despite both states having the 
same envelope for their interference fringes, and despite the 
NOON state being able to reduce the interference period of 
the fringes by a factor of N, the peak multiphoton absorption 
rate of a NOON state is lower than that of a classical state by a 
factor of 2 N ~ 1 . In Sec. [IV] I shall discuss the formalism in the 
paraxial regime, where it is acceptable to regard photons as 
spatial objects described by a configuration-space probability 
density. I shall then investigate the proof-of-concept quantum 
lithography experiment by D' Angelo et al. lfl3ll and show that 
the experiment requires a condition that necessarily reduces 
the two-photon absorption rate, in order to emulate quantum 
lithography accurately. Finally, I shall study the multiphoton 
absorption of a jointly Gaussian multiphoton state, in order 
to illustrate the trade-off between resolution enhancement and 
multiphoton absorption rate. 



The electric field operator is then given by [ 1 

E(x,y,z,t) = E ( - + \x,y,z,t)+h.c 



(3) 



E (+} (x,y,z,t) 



(23T)a sJ V2eo 
x a(k x ,ky,k z ,s)e(k x ,ky,k z ,s) 
x exp(ik x x + iky.y + ik z z — icot), (4) 



where h.c. denotes Hermitian conjugate, and e(k x ,k y ,k z ,s) is 
the unit polarization vector corresponding to one of the two 
polarizations orthogonal to the wave vector. In the following I 
shall consider only the modes that propagate in the positive z 
direction in the z —x plane, and only the s polarization normal 
to the z — x plane. This is consistent with Boto et al. 's pro- 
posal, and equivalent to assuming k z > 0, k y s a and choosing 
one s such that e = y. Following Blow et al. Q, I shall make 
the following substitution to neglect the y dimension: 



dk Y 



2k 

v 



a(k x ,k y ,k z ,s) — > a(k x ,k z )\ 
2(k x ,k z ),a\k' x ,k' z )} = 8(k x -k' x )8(k z -k' z ) 



(5) 

(6) 
(7) 



where L y is the normalization length scale in the y dimension. 
The electric field is thus simplified to 



?{+) 



x,z,t) 



dk.i dk z i—Y 

In^/LyJ-™ Jo \2eoJ 
x a(k x ,k z )exp(ik x x + ik z z — icot). (8) 



II. QUANTIZATION OF TWO-DIMENSIONAL, 
S-POLARIZED, APPROXIMATELY MONOCHROMATIC 
ELECTROMAGNETIC FIELDS 

A. Two Dimensional Approximation 

I shall start with the most general commutation relations for 
creation and annihilation operators of continuous electromag- 
netic field modes in free space [12]: 

[a(k x , k y ,k z ,s) , a 1 " (k' x , k' y , k' z , s')} 
= 8(k x - k' x ) 8 (k y -k y )8(k z - k' z )8 ss , , (1) 

where k x , ky, and k z are the independent continuous variables 
for each mode of the electromagnetic fields, and s denotes one 
of the two polarizations perpendicular to the wave vector. The 
dependent variable in this case is frequency CO, determined by 
the dispersion relation 



B. Propagating Fields 

To conform to classical optics conventions, I shall make k x 
and co the independent variables and k z the dependent vari- 
able, following the procedure of Yuen and Shapiro [15]. This 
coordinate transformation yields 

CO 

dk x dk z = -^—dk x dco, (9) 

c kz 

, i f m l c r ( hco\? co 

(x,t,z = 1=] dk x \ dcoi — ) - f - 

2K^jLy)-mlc Jo \2eoJ C 2 k z 

x a (k x ,k z ) exp(ik x x— icot + ik z z), (10) 

where k z = y/ co 2 /c 2 — k\ is now the dependent variable. Con- 
sider the commutation relation for a(k x , k z ) in terms of the new 
variables, 

[a{k x ,k z )^KK)\ = 8{k x -k' x )8{k z -k' z ) (11) 

= c -J^8{k x -k! x )8{(0- to'), (12) 



co 2 = c 2 {k 2 x +kl + kl). 



(2) 



where the factor c 2 k z /co comes from Eq. (|9]). A new annihila- 
tion operator in terms of k x and co should therefore be defined 
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as [1 



co \ 1 

d(k x ,co) = ( ) d(k x ,k z ), 



(13) 



so that the new operators have the desired commutator, 

[a(k x , co), a f (4, co')] = 8(k x -kf x )8((0- co'). (14) 
Writing k x as K as a shorthand, the electric field becomes 

/ ft \ 2 /-co/c 

£ (+ ^(jc,f,z) = i( — — —I / dK I dco 



8n 2 eoc 2 L y J J- m / c Jo 

CO 



(C0 2 /c 2 -K 2 y/* 

x a (fc, co) exp (r'Kx — i'co/ + ik z z) , (15) 

which is now expressed in terms of z-propagating modes, with 
transverse momentum K and frequency CO as the independent 
degrees of freedom. 

C. Monochromatic Approximation 

We have now obtained a formalism that treats k and CO as 
independent degrees of freedom and corresponds to the exper- 
imental situation of quantum lithography, where optical fields 
are considered as propagating modes. The spatial quantum 
properties of the propagating waves are thus independent of 
the temporal properties. In order to study only the spatial ef- 
fect of resolution enhancement on the multiphoton absorption 
rate, separate from the temporal effects studied by Javanainen 
and Gould and Perina et al. I shall assume, in consis- 
tency with Boto et a/.'s formalism, that the photons are all 
approximately monochromatic with the same frequency CO. 
Again, following the conventions of Blow et al. 1 1411 . 

jdco^Y' a ( K '®)^"( K )Jl^, ( 16 ) 
[d(K) ) a^K')} = 8(K-K'), (17) 

where T is the normalization time scale. The electric field 
envelope is then defined as 



E {+) (x,z, t) = E (+) (x,z) exp{~icot ) , 
i 

i 



(18) 



= i ( ^-V / dKY(K)d(K)exp(iKx + ik,z) , 
\2nJ J-io/c 

(19) 



where 



I? = 



hco 



2eocL y T ' 



(20) 



the magnitude of which is on the order of the optical intensity 
per unit length in x for one photon. y(jc) is a geometric factor, 



y(?c) = 



1 



{\-c 2 K 2 /co 2 yi^ 



(21) 



which arises owing to the invariance of the formalism with 
respect to rotation in the z — x plane. See Appendix lAl for a 
detailed discussion on the physical significance of j(k). 
D. Continuous Fock Space Representation 

With the commutator described by Eq. (fTTT i, and the elec- 
tric field envelope in terms of the photon annihilation oper- 
ator in Eq. ( fT9l , a rigorous quantization of two-dimensional, 
^-polarized, approximately monochromatic optical fields has 
been established. To account for all possible configurations of 
a Fock state in terms of the continuous transverse momentum, 
I shall define the following -photon eigenstate lfl2[[l7ll . 



\K U ...,K N ) 



fm 



? + (jC 1 )...fl t (K- W )|0). 



(22) 



A momentum-space representation of a Fock state \N) can be 
written as IU7I1 



0(Ki, 



,K5Vj 



(jq,...,K5v|A0 
1 



/AH 



(23) 

a(Ki)...a(K N )\N), (24) 



and the Fock state can then be written in terms of this repre- 
sentation, 



rca/c 

\N) = d.K\. 

J -a>/c 

rco/c 

dK\. 

-m/c 



io/c 

-m/c 
co/c 

-m/c 



dK N \K U -,K N )(Ki,...,K N \N) (25) 



dK N (j)(KU:.,K N )\K u ...,K N ). (26) 



(j) is hereby defined as the momentum-space multiphoton am- 
plitude, which describes the configurations of fc's for Af pho- 
tons. <p must evidently satisfy the normalization condition, 



m/c 



dK] 



m/c 



dK N \<j>(Ki )) ...,K N )\ 2 = l, (27) 



' — co/c J-m/c 
and the boson symmetrization condition, 
0(. ..,«•„,. ..,K m ,...) = </>(...., >c m ,...,K-„,...) for any n,m. (28) 

E. Af-Photon Measurements at the Observation Plane 

We shall now observe the photons at z = 0, define (x) = 
E^ + \x,Q), and a spatial multiphoton amplitude \j/(x\ , ...,xn) 
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(0\£(+\x l )...£( + \x N )\N)=VW.i N T 1 Ty(x 1 ,...,x N ), 



yr(xi,... i Xft) = 



(2n) N / 2 J-m/c 



m/c 



m/c 



m/c 



dK N y{K N ) 



X 0(jCi,...,K3v)exp ■ 

The physical significance of y/ is that its magnitude squared is proportional to the ideal Af -photon coincidence rate, 

(:t( Xl ).J(x N ) :) = (N\e(-\ Xi )...e(-\x n )E^( Xi )...eM( Xn )\N) 
= (0\E^( Xl )...E^\x N )\N)\ 2 
= N\t] n \y(xi,...,x n )\ 2 , 



where 



I(x)=E^\x)E^ +) (x) 



(29) 

(30) 

(31) 
(32) 
(33) 

(34) 



is the optical intensity operator. Classically, an ideal A^-photon absorption pattern is given by I (x). In quantum optics, the 
average A^-photon absorption rate becomes 



(:t N ( X ) :) = (:/».../(*) :) = Nlr 1 N \ W (x,...,x)\ 2 , 
which is the Af -photon coincidence rate evaluated at the same position x\ = ... = xtf =x. 



(35) 



F. Upper Bound of A'-Photon Absorption Rate for an W-Photon State 

With the formalism outlined above, it turns out that one can already derive an upper bound for the Af -photon absorption rate 
of an Af -photon Fock state, without knowing the specific form of 0, using Schwarz's inequality | (f\g) \ 2 < (f\f) (g\g), 



\\j/( Xl ,..., XN )\ 2 



1 



m/c rm/c I N 

. dKrjiKx)... I dK N y(K N )<j>(Ki,...,K N )exp[iYK n x n 

(2k) n I 1 J-m/c J-m/c \ ^ 



< 



(2nY 



m/c 
-m/c 



dK\. 



m/c 
-m/c 



dK N \<p{K U ...,K N )\ 7 



m/c r®l c 
dK\... I dx N 

m/c J-m/c 



N 

n 

n=l 



Y( K n) ex P( iK nX„) 



< 



[2nf 



m/c 
-m/c 



dK\y(K) 



7t\N 
I, 



(36) 



(37) 
(38) 



where A — 2nc/(£) is the free-space wavelength. Hence, the 
Af -photon absorption rate has a upper bound, 

(■^{x) :) =Nlr 1 N \ V ( x ,..., x )\ 2 <N\ (^f. (39) 

Recall that tj is on the order of the one-photon optical inten- 
sity per unit length in X . The upper bound shows that the best 
multiphoton absorption rate, regardless of the form of 0, is on 
the order of Iq, where Iq is the optical intensity of one photon 
focused to a width A. Although this upper bound is derived 
for the simple case of two-dimensional monochromatic opti- 



cal fields focused in one dimension, one expects that the situa- 
tion should remain qualitatively similar when the y dimension 
is also considered, leading to a maximum absorption rate on 
the order of the Iq, where Iq becomes the intensity of a photon 
focused to an area of A 2 . The enhancement of the multiphoton 
absorption rate using nonclassical spatial properties of pho- 
tons, if any, is therefore likely to be very limited, compared 
with the linear dependence of the absorption rate on Iq ob- 
tainable using frequency-anticorrelated photons J9I IToh . This 
is due to the resolution limit in the spatial domain that limits 
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the spatial bandwidth of the optical fields, as well as the per- 
fectly local spatial response of Af -photon absorption assumed 
in Eq. $35^. 



III. MULTIPHOTON ABSORPTION RATE OF QUANTUM 
LITHOGRAPHY 



The chief results of Sec. [TT] applicable to quantum lithogra- 
phy are the definition of a normlizable momentum-space mul- 
tiphoton amplitude, Eq. i24\ . which is able to describe arbi- 
trary configurations of quantized, approximately monochro- 
matic, two-dimensional, s-polarized optical fields containing 
N photons, the definition of a spatial multiphoton amplitude, 
Eq. (l30l l. and the average Af -photon absorption rate, Eq. d35l ). 
in terms of the spatial amplitude. In the following I shall use 
these results to calculate the Af -photon absorption rates for a 
NOON state and a classical state with otherwise identical pa- 
rameters. 



A. TV-Photon Absorption of a NOON State 




|/V>|0>+|0>|Af> 





FIG. 1: (Color online). Schematic of quantum lithography by the use 
of a NOON state. 

In its simplest and most essential form, quantum lithog- 
raphy entails the A^-photon absorption of a NOON state [4] 



(Kg. ID, 



\NOON) = -j= (\N) a \0)b + \0)a\N)b) (40) 
= -±=[(Pf + |0), (41) 

where A and B label the two interfering beams (Fig.[T|), and A" 
and & are the creation operators for the two arms. In the con- 
tinuous momentum space, I shall express the corresponding 
annihilation operators as 



, 1 / K+K \ , 

dK-=f\ — — U(k), 



'Ak 



B= dK 



'Ak 



f - 



Ak 

K- Kq 

Ak 



(42) 
(43) 



where f(q) is a normalizable function of a dimensionless pa- 
rameter q, which satisfies / dq\f(q)\ 2 — 1 and describes the 
momentum spread of modes A and B. Ak is the momentum 
bandwidth, and K is the tilt of the two arms. f((K+ k )/Ak) 
and /(— (k — Kn) /Ak) are also assumed to be orthogonal, 



K- Kq 



0. 



(44) 



Ak J \ Ak 
so that A and B satisfy the commutation relations 

\A,A*] = [8,&] = 1, [A,#] = 0. (45) 
The momentum space amplitude, according to Eq. ( l24l i. is 
1 



<1>noon(ki,...,k n ) 



:(0\a(K l )...a(K N )\NOON) (46) 



1 



V2AK iV 

N 

+Uf 

n=\ 



N 

w 



K„ + Kq 
Ak 



Kb 



Ak 



(47) 



One can check that this amplitude satisfies the normalization 
condition, Eq. ( f27b . y/jvocw is thus determined to be 



Wnoon{ x ii—i x n) 



^2(2kAkY 



exp | -z'Kb £ x„ j f\JdK n Y(K„ - Kb)/ exp(/K„jc„) 



n=l I n— 1 



+ exp I iKb £ x » ] EI / dK n7i K n + Kb)/ \-J^J exp(/K„^„) 



(48) 
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At x\ = ... — xn = x, \j/NOON becomes 
Wnoon(x,-.,x) = 



y/2(2KAK) N 

■f exp (iNKqx) 



exp(— iNKox) 



dKj{K- Kb)/ ( — J exp(iK-x) 



dKy(K+ K {) )f (~^) exp(iKx) 



(49) 



r 



I shall define a beam envelope function 

F(x) ee -J= J dK Y (K- Kb)/(^) exp(toc), (50) 

which is simply the electric field envelope of one of the optical 
beams. We then have 

Vnoon(x,- ..,x) = -^=[F N (x)exp(-iNK x) 



■F N (-x)exp(iNKoi 



(51) 



Assuming for simplicity an appropriate /(k/Ak) such that 
F(x) is even, Ynoon is further simplified to 

Vnoon(x,...,x) = \flF n (x)cos(Nk x), (52) 
and the Af -photon absorption rate becomes 

(: I N (x) :) m0N = 2N\n N \F{x)\ 2N c , 2 {NK,x). (53) 

The pattern consists of an envelope \F (x)\ 2N and an interfer- 
ence fringe pattern cos 2 (A^Kb-^)- If AfC << Kq, the envelope 
is much broader than each period of the interference fringes, 
we then obtain the main result derived by Boto et at, which is 
a multiphoton interference pattern cos 2 (Nkqx) with a period 
equal to k/{Nkq) and inversely proportional to N. 

B. A'-Photon Absorption of a Classical State 




FIG. 2: (Color online). Schematic of classical multiphoton lithogra- 
phy. 



The NOON state should be compared with a classical N- 
photon state given by 



11V 



l {A^+ff 



/AH 



V2 



10), 



(54) 



which can be obtained, for example, by putting an Af -photon 
state to one of the inputs of a 50%-50% beam splitter (Fig. [2]). 
The momentum space amplitude is 



N j 

0c(ki,...,K3v) = YI~t= 

b=1 V2 



K„ + Kq 

Ak 



+ f - 



Ak 



(55) 



The amplitude is a product of one-photon amplitudes, under- 
lining its classical nature, y/ becomes 



N J 

\j/ c (xi,...,x N ) = Y[—7=[F(x n )exp(-iK x n ) 

„=i V2 

+ F(-x„)exp(iK x n )]. 



(56) 



Assuming again that F(x) is even, the Af -photon absorption 
rate is 

(■^(x) :) c =2 N N\r 1 N \F(x)\ 2N C o S 2N (K x), (57) 

which has the same envelope ^(x)! 2 ^ as the NOON state, 
although the interference period is fixed at 7z/kq. The peak 
Af -photon absorption rate is 

(58) 



(■J N {0) :) c = 2 n NIt 1 n \F[Q)\ 



2N 



Thus, even though we have meticulously carried out quan- 
tization and normalization, we find that, under very general 
conditions, the peak multiphoton absorption rate of a classical 
state is higher than that of a NOON state by a factor of 2 N ~ 1 , 
despite both having the same envelope ^(jc)! 2 ". Hence, al- 
though the NOON state is able to offer an A^-fold enhance- 
ment of multiphoton interference resolution, the NOON state 
manifestly does not have an enhanced multiphoton absorption 
rate, as claimed by Boto et al. 



IV. PARAXIAL REGIME 

The factor y(fc) makes analytic calculations of the multi- 
photon absorption pattern more difficult, and prevents an in- 
tuitive understanding of the trade-off between resolution en- 
hancement and multiphoton absorption rate. To mitigate this 
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issue, in this section I shall work in the paraxial regime, where 
K <C co/c, and y(>c) 1. This approximation simplifies the 
analysis significantly and adequately describes most optics ex- 
periments, including the proof-of-concept quantum lithogra- 
phy experiment by D' Angelo et al. [13]. 

To justify the paraxial approximation, consider Fig. [3j 
which plots y(fc) with respect to the normalized parameter 
NA = ck/co, defined as the numerical aperture 11811 . One can 
see that y(?c) is relatively flat and w 1 for a wide range of K. 
Even for an NA as high as 0.8, y(jc) is only approximately 1 .3, 
so in most practical cases, y( k) provides only a qualitatively 
unimportant correction factor. 



A. Simple Model of Multiphoton Absorption 



-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

NA = ck/lu 

FIG. 3: (Color online). A plot of j(k) versus the numerical aperture 
NA — ck/co. 

In the paraxial regime, \]f defined in Eq. ( f30b becomes the 
familiar A^-dimensional Fourier transform of 0, 



\jf(xi,...,x N ) 



(2tc) n / 2 



dK\. 



x 0(jCi,...,jQv)exp ( i£ K " x ") ' ( 59 ) 



because k„ <C co/c and j{K n ) 
normalized, 



1 . y/ is then approximately 



/ 



dx l ...dx N \y/(x u ...,x N )\ 2 s» 1. 



(60) 



\\l/(xi, ...,xn)\ 2 can thus be roughly regarded as the 
configuration-space probability density of finding N photons 
near positions x\ , ...,xjy respectively. Provided that we do not 
localize them too precisely, photons as particles in space are 
hence an acceptable concept in the paraxial regime and de- 
scribed by a properly normalized probability density. 

The configuration-space model has been successfully ap- 
plied to the quantum theory of optical solitons fl9ll . where the 
slowly varying temporal envelope approximation holds, so it 
is perhaps not surprising that the model can also be applied 
to the spatial paraxial domain, where the optical beam is rela- 
tively uniform. 




FIG. 4: (Color online). A simple model of a multiphoton absorption 
material, consisting of many small individual multiphoton absorbers. 

This probabilistic spatial interpretation of photons also pro- 
vides an intuitive understanding of the expression for the mul- 
tiphoton absorption rate in Eq. d35l l. Consider a toy model for 
an -photon absorption material consisting of individual N- 
photon absorbers, each occupying a width of , as depicted 
in Fig. 2] The probability of all photons hitting the rath ab- 
sorber situated at £, m , thus exciting an A^-photon absorption 
event at this absorber, is given by 



&»+A£/2 
€m-A§/2 



dx\ 



6n+A§/2 
6n-A?/2 



dx N \^f{x u ...,x N )\ 2 



(61) 



where P(£, m ) is the probability density of the A^-photon ab- 
sorption event. The spatial resolution of the N-photon absorp- 
tion pattern evidently depends on A| . To eliminate this depen- 
dence and make the resolution depend solely on the resolution 
of the optical fields, we shall make AS, very small, 



P^ m )A^A^\ ¥ ^ m ,...,^„ 



(62) 



so that in the limit of a continuous Af -photon absorption mate- 
rial, the probability density of Af -photon absorption becomes 



P(x)=A^- 1 \ W (x,...,x 



(63) 



which is proportional to (: I (x) :) given by Eq.d35b. So, intu- 
itively, an Af -photon absorption event occurs when all photons 
arrive within a very small neighborhood, and the multipho- 
ton absorption pattern is therefore approximately given by the 
conditional probability distribution when all photons arrive at 
the same place. This model has been used by Steuernagel to 
approximate a four-photon absorption material by four dis- 
crete detectors [11], although the explicit derivation here by 
the use of a configuration space model confirms the intuition 
that a multiphoton absorption event occurs when all photons 
arrive within a small neighborhood. 

It must be stressed that although the interpretation of 
\\l/(x\, ...,xn)\ 2 as a configuration-space probability density is 
only valid in the paraxial regime, the expression Eq. (|35| > is 
always a valid description of an ideal multiphoton absorption 
process, because of its dependence on the optical intensity, a 
physically measurable quantity. 
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B. Analysis of Proof-of-Concept Experiment by D'Angelo et al. 

The proof-of-concept quantum lithography experiment by 
D'Angelo et al. II 1311 remains well within the paraxial regime, 
so an explicit analysis of the results can be carried out rela- 
tively easily. In this section, I shall show that the coincidence 
rate detected by in D'Angelo et al. 's experiment is necessarily 
reduced, in order to emulate quantum lithography accurately. 



|0>+£ |2> Far Field 



Pump 



X 



(2) 



1° 

I 



<:f(6):> 



FIG. 5: (Color online). Schematic of D'Angelo et al.'s proof-of- 
concept quantum lithography experiment ITHl . 

In the experiment depicted by Fig. [5] the spontaneously 
generated photon pair has the following quantum state, 

\V)tv\0)+e\2), (64) 

where £ depends on the efficiency of the spontaneous para- 
metric down conversion process, and must remain <C 1, so 
that the quantum state contains only zero or two photons in 
most cases. In D'Angelo et al.'s analysis, the photon pair im- 
mediately exiting the crystal is assumed to have perfect anti- 
correlation in transverse momentum, 



<t>(Ki,K 2 ) = -^=(k u k 2 \2) ~ 8(Kl + K 2 ), 



(65) 



as the pump beam is assumed to be relatively uniform across 
the transverse plane of the crystal and the crystal is relatively 
short. The spatial biphoton amplitude becomes 



Y(xi,x 2 ) ~ 5(*i -x 2 ) 



(66) 



and the photons are assumed to be perfectly correlated in 
space. Clearly, both Eq. d65l l and Eq. ( |66l l are approximations. 
In any case, we shall first follow the approximate analysis and 
normalize the expressions later. Immediately after exiting the 
crystal, the photon pair passes through two slits of width a 
spaced b apart, resulting in a spatial amplitude 



V(xi,x 2 ) ~ 8{x\-x 2 ) 



2 

n 



rect 



■b/2 



-rect 



-b/2 



(67) 



8(x\-x 2 ) 



rect 



X\ +x 2 — b 
2a 



rect 



x\ +x 2 - 
2a 



(68) 



Because the photons are assumed to be perfectly correlated 
in space, they always pass through the same slit, resulting in 



a NOON state in the spatial domain. The momentum-space 
amplitude, on the other hand, is given by 



0(ki,K2) ~ sine 



a(Ki + K 2 ) 



cos 



b(K\ + K 2 ) 



(69) 



This is obviously not the NOON state in the momentum space 
for quantum lithography. To emulate quantum lithography in- 
directly, D'Angelo et al. then let the photons propagate to the 
far field. Via Fraunhofer diffraction lfl8ll . the angular two- 
photon coincidence distribution, (:P(d) :) 112011 . becomes the 
magnitude squared of the Fourier transform of \j/(xi,x 2 ), or 



<:/ 2 (0):>oc| e |2 



x [2nd 2n9\ 



,, 2 . 2 [2Kad\ , [2nbd\ 
\e\ smc f—jco. (^-j 



(70) 
(71) 



This expression is the same as that derived in Ref. IIOl . |e| 2 is 
now regarded as the total probability of two photons reaching 
the detection plane. To be more rigorous, however, \\r(x\,x 2 ) 
in Eq. d68l > and <$>{K\,K 2 ) in Eq. d69l ) need to be normalized. 
For example, the delta function in Eq. d68l ) should be replaced 
by a sharp normalizable function, 



8{x\ 



\ 1 / x\ 



x 2 



a 



(72) 



where g(q) is a function of a dimensionless parameter q and 
is normalized according to / dq\g(q)\ 2 = 1. a is defined as 
the biphoton coherence length, which depends on the phase 
matching condition of the parametric down conversion pro- 
cess and the nonlinear crystal length, a is assumed to be much 
smaller than a and b, but must still be non-zero in reality. The 
normalized y/ then becomes 




rect 



and the normalized (j) becomes 

^oca _ \a(Ki - k 2 ) 



x\ +x 2 + b 
2a 



(73) 



lj>[Ki,K 2 ) = 



a(K\ + K 2 ) 



b{Ki + K 2 ) 



where 



G(P) 



dqg(q)exp(-ipq) 



(74) 



(75) 



is the dimensionless Fourier transform of g. The angular dis- 
tribution is hence 

(■ I (8) ■) - |e| ^2 |G(0)rsinc- ( J cos 2 



J' 
(76) 
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which is proportional to a. Thus, in order to produce a NOON 
state in the near field and emulate quantum lithography accu- 
rately in the far field, the photons need to pass through the 
same slit, the biphoton coherence length needs to be small 
and is even assumed to be zero in the analysis by D' Angelo et 
ai, but then the coincidence rate, proportional to the biphoton 
coherence length, is necessarily reduced. 



C. Jointly Gaussian Multiphoton State 

In Sec. |nl] we have studied the use of NOON state for 
quantum lithography, and it has been shown that the NOON 
state has a lower multiphoton absorption rate than a classi- 
cal state. In Sec. IIV Bl we have also seen that, in order to 
approximate a NOON state accurately in D' Angelo et al.'s 
experiment, the multiphoton absorption rate is necessarily re- 
duced. While these results provide evidence that it is proba- 
bly impractical to use a NOON state for multiphoton lithogra- 
phy, the NOON state is only one example of infinitely many 
possible quantum states for optical fields, and other quantum 
states might be able to perform better while still producing 
an enhanced resolution. For example, instead of producing 
enhanced interference fringes with a minimum period on the 
order of X/N, Bjork et al. 112 111 considered another special 
quantum state, called the reciprocal binomial state, in order 
to produce a sharp interference spot, with a minimum width 
on the order of A /N, within a periodic pattern. Still, it remains 
a question whether this state can produce a significantly better 
multiphoton absorption rate, as the NOON state is still a sig- 
nificant component of the reciprocal binomial state. Steuer- 
nagel, in particular, studied the four-photon reciprocal bino- 
mial state, with four discrete detectors approximating an ideal 
four-photon absorption material, and found that the multipho- 
ton absorption rate is worse than that of a classical state lUlll . 

In this section, I shall study an arguably simpler and more 
intuitive iV-photon state that produces a quantum-enhanced 
Gaussian multiphoton absorption spot, in the paraxial regime. 
I shall call this state a jointly Gaussian state, which is a 
quantum generalization of the well known classical Gaussian 
beams and is able to account for quantum correlations of the 
photons. It is shown that, in certain limits, the jointly Gaus- 
sian state is also able to reduce the size of the multiphoton 
absorption spot by a factor of N compared with the one pho- 
ton case, but the reduction of size is always accompanied by 
a reduced multiphoton absorption rate. For the jointly Gaus- 
sian state, the quantum correlations of the photons, the size 
of the multiphoton absorption spot, as well as the absorption 
rate can all be adjusted by changing just two parameters, so a 
study of this state is able to quantify and elucidate the trade-off 
between resolution enhancement and multiphoton absorption 
rate. 



used in many-body physics, which will significantly simplify 
the analysis later. This coordinate transformation is 



1 N 

K = Tp T" 1 , Ki: K ' n = K n~ K; 



(77) 



where K is the average momentum, and K' n is relative momen- 
tum. K and N — 1 of the K' n 's form a complete basis, so I shall 
somewhat arbitrarily define k' n as the extraneous linearly de- 
pendent variable, 



The new differential is 



N-l 



1 



dKdK , l ...dK , N _ l = —dK\...dK N , 
so a new multiphoton amplitude should be defined as 

<I>'(K,k / 1 ,...,k' n _ 1 ) = Vn<I>(K + k[,...,K+k' n ), 
and normalized as 



dKdK , 1 ...dK , N _ 1 \<j)'(K,K' 1 ,...,4f_ i 



= 1. 



(78) 



(79) 



(80) 



(81) 



The multiphoton absorption rate in terms of the new amplitude 
in the paraxial approximation is 



y/N J dKdK , l ...dK , N _ l 



x <j>'(K, k(, k' n _ , ) exp (iNKx) 



(82) 



The multiphoton absorption rate is therefore the magnitude 
squared of the one-dimensional Fourier transform of <$>', with 
respect to only the total momentum NK. 



2. N -Photon Absorption of a Jointly Gaussian State 

In terms of the new coordinate system, we can now define 
the jointly Gaussian state as follows, 



f (K,k[,...,k' n _ ] 



vCexp I 



AB 2 



exp 




(83) 



where B and j3 are two parameters assumed to be real for sim- 
plicity, C is the normalization constant, 



1. Many-Body Coordinate System 

Before defining a jointly Gaussian state, I shall first take a 
detour and define a new many -body coordinate system, widely 



C = 



N 



(2nY 



1 



BP 



N-l 



(84) 



as derived in Appendix IB II This definition is inspired by 
the well known jointly Gaussian distribution in statistics l22tl . 
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The form of Eq. (1831 is much simpler than a general jointly 
Gaussian distribution because of bosonic symmetry, as dis- 
cussed in Appendix IB 21 The momentum amplitude of the 
state given by Eq. d83l is also very close to that of a soliton 
state 11911 . so one can obtain a general jointly Gaussian state 
approximately by adiabatic control of spatial solitons |23]. 
The covariances of K and are calculated in Appendix 



B 



(85) 



so B is a measure of the spread in the average momentum, and 
j3 is a measure of the spread in the momentum relative to the 
average. 

The variance of K„, the momentum of each photon in the 
original coordinates, is also derived in Appendix IB 41 and 
given by 



(86) 



K n must be smaller than co/c, otherwise k z would become 
imaginary, and K n must be much smaller than co/c for the 
paraxial approximation to hold. Moreover, if an optical sys- 
tem has a certain aperture, it would also limit the transverse 
spatial frequency 111 811 . Hence the variance of K n , (^ 2 ), must 
be limited, and there exists a trade-off between B and ft . 

The multiphoton absorption pattern of the jointly Gaussian 
state can be determined using Eq. 



(:r(x):) = N\T 1 N VN 



B/T^expl 



-2N 2 B 2 x 2 ) 



(87) 



The pattern is a Gaussian, with a root-mean-square width 
given by 



W 



Jdxx 2 (:I N (x)-.y 
Jdx(:I N (x) :) 



1 



4NB 



(88) 



First, consider the case in which the photons are uncorrected, 
and the classical Gaussian state in the original system of coor- 
dinates is given by a product of one-photon Gaussian ampli- 
tudes, 



0c(jCi,...,JQv)' 

As shown in Appendix 
Gaussian state when 



A' 

ri ex p 

«=i 



(89) 



this corresponds to the jointly 



7 P 2 

B — — 



32 
N 



N 



(90) 



The classical variance of the average momentum is equal to 
the variance of each momentum, (k 2 ), divided by N. This 
is consistent with the statistics of independent photons. The 
multiphoton absorption width becomes 



W C = 



1 



(91) 



where the subscript C denotes the value for a classical state. 
Equation (l9TT l can be regarded as the standard quantum limit, 
and is better than the one-photon case by a factor of \/N. On 
the other hand, the mi nimum width is obtained when we max- 
imize B so that B = y/ (k%) and let /3 be zero, 

1 



Wm 



(92) 



The factor-of-A^ enhancement compared with one-photon ab- 
sorption, or the factor-of-VAf enhancement compared with 
classical Af -photon absorption, can be regarded as the ultimate 
quantum limit, and is consistent with other quantum enhance- 
ment schemes j24Tl . This enhancement, however, comes with 
a heavy price. Equation d87l l shows that the multiphoton ab- 
sorption rate is proportional to Bfi N ~\ so while increasing B 
and reducing j3 makes the Gaussian pattern sharper, the reduc- 
tion in j3 also reduces the multiphoton absorption rate, more 
so for large N. 

To quantify this trade-off, we shall fix (k 2 ") as a given re- 
source, and define a spot size reduction factor r with respect 
to the classical case, 



(93) 



so that r = 1 corresponds to the standard quantum limit, and 
r = \/N corresponds to the ultimate quantum limit. We shall 
also define a normalized peak absorption rate R with respect 
to the rate in the classical case, 




R 



_ (: ^(0):) 
<:/*(0):> c 



N-r 1 
N- 1 



(94) 




— N = 


2 


---N = 


4 


— N = 


9 


■-N - 


00 
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Spot Size Reduction Factor, r 

FIG. 6: (Color online). Plots of peak multiphoton absorption rate 
versus spot size reduction for several values of N. Both quantities 
are normalized with respect to classical values. Interestingly, in the 
limit of N -> 00, R -> r exp [ ( 1 - r 2 ) /2] . 

Figure |6]plots R versus r for several values of N. This result 
is decidedly disappointing, as it shows that the maximum mul- 
tiphoton absorption rate is obtained when the state is a classi- 
cal state, or r = 1, and the peak rate monotonically decreases 
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to zero as the spot size is reduced. Furthermore, even if one 
is willing to sacrifice the resolution and increase j3, the peak 
absorption rate is still reduced, because of its dependence on 
B. 

For applications such as multiphoton spectroscopy, spatial 
resolution is not important, and it is more desirable to maxi- 
mize the total multiphoton absorption rate. We can define the 
normalized total rate as 



Rt, 



Jdx(:I N {x) :) 
/<&<:/*(*) :) c 



N- 



N-l 



(95) 



N = 


2 


N = 


4 


N = 


9 
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FIG. 7: (Color online). Plots of total multiphoton absorption rate 
versus spot size reduction for several values of N. 

Figure|7]plots R tot versus r. It can be seen that the total rate 
does increase when one increases the spot size, but the rate 
enhancement is very moderate. In fact, in the limit of N — > °°, 
Rtot approaches 



Rt, 



exp 



1 



(96) 



so the ultimate rate enhancement, when resolution is com- 
pletely sacrificed and r = 0, asymptotically approaches 
exp(0.5) w 1 .65 for large N. This small enhancement of mul- 
tiphoton absorption rate is not likely to be useful. 

To understand the above results, it is helpful to consider the 
position correlation of the photons, derived in Appendix[ 



(x n x n 



-(—-- 
4N \NB 2 P 2 



n/m. 



(97) 



To obtain an enhanced multiphoton resolution, the bandwidth 
of the average momentum, B, must be increased, leading to a 
positive correlation in the momenta and a negative correlation 
in the positions of the photons. So the photons are actually 
less likely to arrive near one another, leading to a lower prob- 
ability of the photons hitting the same absorber and therefore 
a correspondingly less multiphoton absorption rate. On the 
other hand, if B is reduced, the position correlation becomes 
positive, and the photons are more likely to arrive close to one 
another, leading to a slightly enhanced total absorption rate. 



Ultimately, how close the photons can arrive with respect to 
one another is still restricted by the resolution limit of the op- 
tical system, so the rate enhancement is not significant. 

It must be stressed again that the preceding argument, as 
well as Boto et al.'s heuristic argument about photons con- 
strained to arrive at the same place, are only applicable to the 
paraxial regime, where the positions of photons are relatively 
well defined quantities. The example of the jointly Gaussian 
state shows that, even in the paraxial regime, Boto et al.'s 
heuristic argument is not correct, and the photons are actually 
less likely to "arrive at the same place" when the resolution is 
enhanced. Although this has been shown by Steuernagel for 
the specific example of four-photon reciprocal binomial state 
ifTTll . the study of joint Gaussian state here confirms this fact 
for an arbitrary number of photons. 



V. DISCUSSION AND CONCLUSION 

In summary, through a rigorous study of quantum lithog- 
raphy and the NOON state, an investigation of the proof-of- 
concept experiment by D'Angelo et q/. lll3ll . and an analysis 
of a jointly Gaussian state, I have been unable to find any evi- 
dence, as far as the spatial domain is concerned, that supports 
the heuristic claim by Boto et al. ||4(], namely that the photons 
would be "constrained to arrive at the same place" and the 
multiphoton absorption rate would be enhanced due to spatial 
effects. On the contrary, all examples show that the multipho- 
ton absorption rate is actually reduced, more so for larger N, 
when a quantum state is used to enhance the resolution. 

Admittedly, there are several assumptions involved in the 
analysis, the negligence of time domain effects in particu- 
lar. As Javanainen and Gould [9] and Perina et al. llOll have 
shown, frequency entanglement can enhance the multiphoton 
absorption rate, but this enhancement is likely to be indepen- 
dent of the detrimental spatial effect, which seems to be an 
unavoidable penalty incurred by the resolution enhancement 
effect itself. That said, it remains to be proved whether tak- 
ing time domain into account would be able to eliminate or 
reverse the detrimental spatial effect. 

Moreover, only the NOON state and the jointly Gaussian 
state have been studied in this paper, but the possibility of 
other exotic quantum states being able to enhance the resolu- 
tion while maintaining a respectable multiphoton absorption 
rate cannot be ruled out. In fact, alternative strategies have 
already been proposed to solve the low exposure problem of 
quantum lithography. For example, Agarwal et al. have shown 
that strong nonclassical beams from a parametric amplifier 
can also produce enhanced two-photon interference fringes 
[5], albeit with a background worse than a classical multiple 
exposure technique 01 . Hemmer et al. also proposed the use 
of a narrowband multiphoton absorption material and classi- 
cal light D7D, similar to Yablonovitch and Vrijen's proposal 
Q3D , to achieve the same resolution enhancement as quantum 
lithography. 

In conclusion, in light of the results set forth, the original 
quantum lithography scheme is unlikely to be practical in the 
near future. Nonetheless, it has inspired many ongoing re- 
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search efforts on the elusive goal of beating the optical reso- 
lution limit, and should therefore remain an interesting theo- 
retical concept. 
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APPENDIX A: PHYSICAL SIGNIFICANCE OF THE 
GEOMETRIC FACTOR j(k) 

To understand why the factor y(jc) arises in Eq. (f30b of 
the formalism, consider the electric field envelope given by 
Eq. <QJ>, 

r<»/c 

E [+, (x,z)°< / dKY(K)d(K)exp(iKx + ik z z) , (Al) 

J-co/c 

where k z is the dependent variable given by k z = 
y/ (Q 2 /c 2 — K 2 . We shall leave the form of y(jc) unspecified 
and derive it purely from the fact that the formalism is invari- 
ant under a rotation in the z — x plane. Imagine that the electric 
field profile is rotated anticlockwise in the z — x plane, where 
z is the horizontal axis and x is the vertical axis, by an angle 
9. This is equivalent to defining new coordinates as follows, 

x' =xcos9 +zsin0, z' = — xsin 6 +zcos 9. (A2) 

The envelope becomes 

E^\x',z') 

m/c 

dKy(K)a(K) 

m/c 

exp [iK (cos 9x' + sin 9z) + ik z (— sin 9x + cos 9z') ] 

(A3) 

m/c 

dKj(K)d(K) 

-m/c 

x exp [/(jfcos 9 — k z sm9)x' + /(fcsin0 +k z cos9)z / ] ■ 

(A4) 

If we define the momenta in the new coordinate system to be 

k' = K°cos 9 — k z sin 9 , k' z — ifsinO +£ z cos 9, (A5) 

evidently the new definitions still satisfy the dispersion rela- 
tion 



K ,2 + k? = K Z - + kt = 



a>- 



-2 ' 



More crucially, the coordinate transformation yields 



dK= -^dK', 



(A6) 



(A7) 



so that the electric field envelope becomes 

rC ° /C j.J k z, 
I -m/c 



, \ r m / c k 

£( + V,z')- / dK'-fY(K)a{K)exp(iK'x' + ik' z z!). 

J-mlc k, 



(A8) 



If the electric field is invariant to such a rotation, we should 
be able to define a new momentum-space operator a(jc') such 
that 

, ^ r m / c 
E(+\x',z')°< 

J -m/c 

with the commutator 

[a'{K') 7 d'\K")] = 8{K'-K"). 

Comparing Eq. ( IA8b and Eq. JA9K we have 



t/K'y(K')a'(K-')exp (/r-V + , (A9) 



(A10) 



(All) 



k' z y(fc' 

but the coordinate transformation also restricts the relation be- 
tween fl(fc) and a'(K'), 



«MH^)W(|)W>. 



Combining Eq. (IA1U and Eq. ( IA12I ) yields 
y(jc) 



Y(k') 



(A12) 



(A13) 



For Eq. dA13l ) to hold for any rotation, j{k) must depend on 
K only according to the following, 

r(«)=ckP = {a)2/c2 C _ K2y/4 , (A14) 

where C is an arbitrary constant. Equation dA14b is identical 
to Eq. d2Tb with C = \J co/c. Hence, the factor y(ff) arises 
purely due to the invariance of the formalism with respect to 
rotation in the z~ x plane. Furthermore, the general trans- 
formation rule for a(fc) with respect to a rotation is given by 
Eq. dAllll . or 



fl(fC) 



(A15) 



Ksin9 +k z cos9\ 1 , . 

' a (kcos9 — k z sm9). (A16) 



APPENDIX B: PROPERTIES OF THE JOINTLY GAUSSIAN 
STATE 

In this section I shall derive several properties of the jointly 
gaussian state, 

= ^ ex P ex P ("4^5 t K ") ' (B1) 

where C is the normalization constant, B and /3 are real pa- 
rameters, and k' n is given by — YIIzI K 'n- 
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1. Normalization 

To calculate C, consider the normalization 
C J dKd k[... J exp ("|^) exp ( -^jL £ ^ = 1 

(b: 

CV^fi J dK[...dK' N _^x V (-^p £ icfj = 1 



(B3) 



As K' N — Y*n=l K n' 



N N-l /N-l 



N-1N-1 



n=l n=\ \n—l 



n=l m= 1 



where A nm = 8 nm + 1 and \A nm \=N. Using the normalization 
of a jointly Gaussian function 022I1 . 



1 



JV— 1JV— 1 



dK[...dK l N _ l exp [ E E K-vAnm^m 



^P 2 „=i m= i 



(^) A '- 1 j3 iV - 1 |A„ m |-2 



N-l o N-l 



N 



(B5) 
(B6) 



Substituting Eq. (|B6]i into Eq. (B3]i gives Eq. ( f84l . 



2. Bosonic Symmetry 

I shall now show that Eq. ( IBU is a consequence of enforc- 
ing bosonic symmetry on a general jointly Gaussian function. 
In the original coordinate system, can be determined from 
Eq. 



<j)(K h ...,K N ) 



exp 



4B 2 N ^ 



1 V f _ 1 V 



which can be rewritten as 

0(K- 1 ,...,Kj V ): 



ex P ( ^ E ^"^nm Kr, 



(B7) 



(B8) 



Equation ( IB8t is a general jointly Gaussian function, but be- 
cause of boson symmetry of (j) as prescribed by Eq. ( f28b , B ; ,„, 
must have identical on-axis components, as well as identical 
off-axis components. With some algebra, B nm can be deter- 
mined from Eq. (1B7I >. 



Bnn — 
B nm 



1 



N 2 B 2 
1 



1- 



1 \ 1 



NJ ]8 



2 » 



JV 2 5 2 A?j3 2 



« 7^ m. 



(B9) 
(BIO) 



Any B„„, with identical on-axis components and identical 
off-axis components can be specified using B and /3, so 
Eq. (ED can specify any general jointly Gaussian functions 
with bosonic symmetry. 



3. Covariances 

The covariances of momentum variables should be deter- 
mined from the probability distribution 



K 2 \ 



1 



N-lN-l 



: CeXP ^"2^2 J 6XP ^~2^2 ^ E K '" AnmK ™j ' (B11) 

where A„,„ = 5„„, + 1 is defined in Appendix lB II The variance 
of K is simply given by 



(K 2 )=B 2 , 



(B12) 



while the covariance matrix for K^'s is given by 

(k'X,) = p 2 A m l (B13) 



Because A nm only has two parameters, its inverse, defined as 
C nm = Aj, is easy to calculate and is given by 



C„„ = l-^; C nm = -—, n^m. (B14) 



We thus obtain the covariances, 



N 



, n^m. (B15) 



4. Classical Gaussian State 

The covariances of k„ in the original coordinate system are 
(k 2 ) = ((K+k>) 2 ) = (K 2 ) + (k> 2 ) (B16) 

(B17) 

P 2 



B 2 +[l-^)P 2 , 



(K„K m ) = {K 2 ) + {K'„K' m )=B 2 -^. (B18) 

So the photons are uncorrected when B 2 = fi 2 /N, and 
<p ( K"i , . . . , Kjy) in Eq. (IB8b can be written as 



C N ( K 

C ( J c 1) ..., K5 v)=y^n ex p^-i|2 



(B19) 



a product of one-photon Gaussian amplitudes, and therefore a 
classical state. 
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5. Configuration-Space Multiphoton Amplitude 

In the paraxial regime, the configuration-space multiphoton 
amplitude can be obtained by Fourier transform of Eq. 



yr(*i,...,Xtf) oe exp 



- l x 



(B20) 



and the covariance matrix for the photon positions is (x n x„ 
B m jj4. The position variance is then 



N 2 B 2 



1- 



1 \ 1 



NJ P 7 



(B22) 



which is determined using the well known characteristic func- 
tion of a jointly Gaussian distribution ll22ll . The configuration- 
space probability density is thus 



|y/(>i,...,jQv)| 2 oc exp -2^x„B„jx m J , 



(B21) 



and the covariance is 



= _L(_L__L 

4-N \NB 2 j5 2 



n =/= m. 



(B23) 
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